Noname manuscript No. 

(will be inserted by the editor) 



The Quasi-Steady State Assumption in an Enzymatically Open System 

Ed Reznik 1 3 , Daniel Segre 1 ' 2 ' 3 , William Erik Sherwood 3 4 * 



January 13, 2013 

Abstract The quasi-steady state assumption (QSSA) forms the basis for rigorous mathematical justification of the 
Michaelis-Menten formalism commonly used in modeling a broad range of intracellular phenomena. A critical 
supposition of QSSA-based analyses is that the underlying biochemical reaction is enzymatically "closed," so 
that free enzyme is neither added to nor removed from the reaction over the relevant time period. Yet there are 
multiple circumstances in living cells under which this assumption may not hold, e.g. during translation of genetic 
elements or metabolic regulatory events. Here we consider a modified version of the most basic enzyme-catalyzed 
reaction which incorporates enzyme input and removal. We extend the QSSA to this enzymatically "open" system, 
computing inner approximations to its dynamics, and we compare the behavior of the full open system, our 
approximations, and the closed system under broad range of kinetic parameters. We also derive conditions under 
which our new approximations are provably valid; numerical simulations demonstrate that our approximations 
remain quite accurate even when these conditions are not satisfied. Finally, we investigate the possibility of damped 
oscillatory behavior in the enzymatically open reaction. 

1 Introduction 

The Michaelis-Menten formalism (TJE) is perhaps the most commonly used framework for modeling the dynamics 
of biochemical processes. Derivation of the formalism proceeds from the assumption that in an enzymatic reaction, 
the level of enzyme bound with substrate to form enzyme-substrate complex rapidly equilibrates at the start of the 
reaction. Following this initial transient, the levels of free and bound enzyme maintain near-equilibrium, or "quasi- 
steady- state," for the remainder of the reaction, during which time their concentrations may be accurately modeled 
as simple rational functions of the substrate concentration. 

The quasi-steady-state assumption (QSSA) has been a topic of mathematical interest for nearly a century 
l2HT2ll. Its modern justification is based upon singular perturbation analyses that exploit separations between the 
time-scales on which reaction variables evolve in order to derive analytic approximations to the reaction dynam- 
ics and conditions under which these approximations are valid [3, 5|. Classical QSSA analyses, like the original 
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Michaelis-Menten derivation, consider biochemical systems which are "closed." That is, they admit neither addi- 
tion nor removal of unbound enzyme or substrate from the reaction once it is underway. This leads to a conserva- 
tion relation which greatly simplifies the mathematical analysis. While this scenario may reasonably approximate 
laboratory experiments, biochemical isolation is not so common in vivo, where neither enzyme nor substrate lev- 
els remain constant. Recent studies have extended the QSSA for enzyme-catalyzed reactions to which substrate 
is added at either at a constant rate (5) or as a periodic input lITQl . The question of the effect of variable enzyme 
levels on the QSSA has remained unaddressed up to now, though enzyme concentrations in living cells continually 
fluctuate due to multiple mechanisms: Transcription and translation of genetic elements may add to the available 
store of enzyme, while degradation (by proteases) and deactivation (for example, by kinases via phosphorylation) 
may subtract from it. Similarly, allosteric regulation by small-molecule metabolites can substantially shrink or ex- 
pand the number of enzyme molecules available and able to catalyze a reaction. In fact, the time-scales on which 
such alterations in available enzyme levels take place may vary widely (from a few microseconds for changes in 
enzyme activation state to minutes for eukaryotic enzyme transcription and translation [[13] [HI). Under a variety 
of different circumstances they may overlap with the intrinsic time scales of the biochemical reactions in question: 
the values for enzymatic reaction rate constants implied by estimates of turnover and reaction velocity span at least 
six orders of magnitude, ranging from microseconds to nearly a second, and these rates themselves are affected 
significantly by pH, temperature, and other factors |[T5l . 

The widespread application of the QSSA and the Michaelis-Menten formalism to the kinetic modeling of 
metabolic events rests in no small part on the common understanding that changes in enzymatic levels occur much 
more slowly than metabolic effects. Relatively minor modifications to the classical QSSA approach may be made 
to address some of the above mentioned scenarios in which the time-scales for changes in free enzyme levels 
and other reaction rates clearly overlap [ 16 1, and indeed, the general success and flexibility of the classical QSSA 
approach has understandably led certain of its premises to remain essentially unexamined. Yet a quantitative 
assessment of the assumption of unchanging enzyme levels is of theoretical and practical interest. How does 
variation in the amount of available enzyme affect the dynamics of biochemical reactions? Can the assumption of 
conservation of total enzyme be relaxed? 

Here we consider an enzymatically "open" biochemical reaction in which both enzyme input and removal are 
allowed. Our goals are threefold: First, we seek to extend the QSSA to this more biologically accurate situation. 
Second, we wish to assess the comparative accuracy with which the dynamics of the enzymatically open system 
are approximated by the enzymatically "closed" system on which classical QSSA analyses are based. Third, using 
the approximations derived in our extension of the QSSA, we aim to describe the range of possible dynamics the 
open system may exhibit in various parameter regimes. 

The remainder of this section briefly reviews the classical QSSA framework and its variants for the closed 
system, drawing primarily upon the exposition in [5 |, and it introduces the mathematical description of the open 
system. 



1.1 The QSSA framework for enzymatically closed systems 

The development of the Michaelis-Menten formalism, as well as the classical QSSA framework for its analysis, 
proceeds by first considering the "canonical" simplest enzyme-catalyzed biochemical reaction, to which we refer 
herein as the closed system or closed reaction: 

S + E^C^>E + P (1) 

k-i 

Substrate S and free enzyme E combine reversibly to form complex SE = C, which irreversibly yields free enzyme 
and product P upon dissociation. Applying the law of mass action, the dynamics of species concentrations for this 
reaction are described by a set of coupled nonlinear ODEs: 
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dE 

— = -k 1 ES+(k- 1 +k 2 )C (2a) 
at 

^ =kiES-(k- X +k 2 )C (2b) 
at 

^- =k-iC-hES (2c) 
at 

dP 

— = k 2 C (2d) 
dt 

The concentration of complex and product is normally taken to be zero at the beginning of the reaction, i.e. 
(E,C,S,P) = (£b,0,5 0j 0) at time t = 0. 

Two conservation relations obtain for the enzymatically closed system: 



S = S(t)+C(t)+P(t) (3a) 
E = E(t)+C(t) (3b) 



Both relations simply express the principle of conservation of mass. Equation ( [3a] ) asserts that substrate is neither 
added to nor removed from the reaction once it has begun (except by conversion to product); Equation ( [3b] ) is the 
analogous assertion for free enzyme, i.e., that the reaction is enzymatically closed. 

The first key step of any classical QSSA-based analysis of ^ is to invoke ( [3b] ) in order to reduce its dimen- 
sionality. Omitting the product formation equation ([2d), which decouples, the reduced system is 



^ =k x (E -C)S-(k-x+k 2 )C (4a) 
at 

^ = k-iC-k x (Eo-C)S (4b) 

The goal of the analysis is then to derive approximate analytic expressions for the concentrations of complex and 
substrate over time. 

The second essential step of the analysis is the recognition that the system Q possesses two distinct time- 
scales, so that its evolution proceeds in two phases. In the first phase, the concentrations of free enzyme and 
complex rapidly equilibrate, while the level of unbound substrate is assumed to remain nearly constant, i.e. S^So. 
During this initial transient (or pre- steady- state) period, the concentration of complex exponentially approaches 
its quasi-equilibrium value 

where S is taken to be So and Km = (k- 1 + k 2 ) jk\ is the Michaelis-Menten constant. Note that the quasi-equilibrium 
is a function of the substrate concentration; this is the basic functional form for Michaelis-Menten dynamics. To- 
gether with conservation relation pa} , Equation §5§ also dictates the quasi-equilibrium level of free enzyme. 

In the second phase, the quasi- steady- state period, the concentration of substrate falls as it is slowly converted 
to product according to the ODE 

dS = k 2 E S 

dt K M + S 1 } 
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The duration of each phase is determined by the characteristic time scale for the species concentration undergoing 
the most rapid change within it. The time scale of the first phase is given by the time constant for the exponential 
relaxation of the complex concentration, holding substrate level constant, 

(7) 



k x (S +K M ) 

A useful estimate for the time scale of the second phase is given by the ratio of the maximum change in substrate 
concentration to the maximum rate of change in substrate concentration, which yields 

ts=^^ (8) 

When the time scales for the two phases are well separated, i.e. t c <O s , and the change in substrate concentra- 
tion during the initial transient period is small, the QSSA is valid, as are consequent analyses which invoke it. 
In particular, the accuracy of singular perturbation approximations to the transient ("inner solutions" or "inner 
approximations") and quasi-steady-state dynamics ("outer solutions" or "outer approximations") of the system 
is assured. The requisite separation of time scales was originally shown to hold when substrate concentration 
greatly exceeds free enzyme concentration, a condition which typically obtains for in vitro biochemical experi- 
ments QJ[3]|5). Subsequent analysis showed that the approximations which follow from the QSSA also obtain 
under opposite circumstances, when the level of free enzyme greatly exceeds that of substrate (5] [7] [IT). QSSA 
which presumes So ^> Eq is usually called standard QSSA, and QSSA which presumes Eo ^> So is typically called 
reverse QSSA. In Section[4]of this paper, we investigate the behavior of an enzymatically open biochemical system 
in both the standard and reverse QSSA regimes. 



1.2 The open reaction: Enzyme input and removal 

The biochemical system studied in this paper is a modification of the canonical closed reaction {T} which admits 
the addition and removal of free enzyme from the system while the reaction is underway: 

0f^E + Sf^C4e + P (9) 

k-3 fc-i 

Free enzyme is introduced de novo to the closed reaction at a constant rate £3 and removed or inactivated (but not 
bound to substrate) at a concentration-dependent rate with rate constant k-3. This we refer to as the open system 
or open reaction. Its dynamics are described by ODEs which differ only slightly from ^ 



dE 
dt 
dC 
~dt 
dS 
dt 
dP 



= k 3 -k- 3 E-kiES+(k- 1 +k 2 )C (10a) 

= k 1 ES-(k- 1 +k 2 )C (10b) 

= k-iC-k x ES (10c) 

= k 2 C (lOd) 



The initial conditions are again taken to be (E,C,S,P) = (Eq,0,Sq,0). 

One possible scenario described by ([9} and {To| ) is a simple enzyme-catalyzed reaction during which new 
enzyme is introduced due to transcription and translation of genetic elements, and also degraded in the cytosol. 
In this case, we expect k 3 and k-3 to be several orders of magnitude smaller than the other rate constants for 
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the reactions. For example, enzyme synthesis takes on the order of minutes in E. coli, while the time needed for 
diffusion-limited equilibrium binding of a small molecules to a protein, which corresponds to k\ in is on the 
order of microseconds to milliseconds |[T3l . However, as mentioned in the Introduction, the range of rate constants 
for enzymatic reactions is broad enough that for many biologically plausible circumstances reasonable values for 
k\,k-\, and &2 may approach the same orders of magnitude as k 3 and/or k-^ (l5\ . 

Though the previous biological example is worth bearing in mind, we emphasize that the biochemical ([9]) 
and mathematical ( [To| descriptions of the open system are generic and might usefully be applied to a range of 
biological situations. The rate of enzyme input, k 3 , may be taken to stand for the rate of increase in availability 
for biochemical participation of free, active enzyme by any of several means ranging from the level of allosteric 
effects to that of gene activity. And similarly, the interpretation of the rate of enzyme removal, k-3, is not strictly 
tied to any specific biological mechanism. Furthermore, our analysis and numerical investigations in the following 
sections do not depend on £3 or k-3 being substantially smaller than other rate constants in the system. Even if one 
discounts the possibility that the time-scales of the enzyme-substrate reaction and enzyme input/removal might 
intermingle in vivo, the framework of ^ and fTQ| ) may be considered as a means by which to examine in some 
generality the QSSA presumption that during biochemical reactions total enzyme levels remain essentially static, 
and in particular, its essentiality for QSSA-based analyses and the consequences for the QSSA were it not to hold. 

The remainder of this paper is devoted to extending standard QSSA analysis to the open system, with particular 
emphasis on dynamics during the initial transient period. In Section [2] we apply the QSSA to derive "inner" 
approximations to the enzyme and complex dynamics of the open system. In Section [3] we calculate time-scales 
for the initial transient and quasi-steady-state dynamics of the open system, and we derive conditions on the 
separation of these time-scales under which the approximations obtained in Section[2]are valid. Section|4]presents 
the results of numerical simulations which we use to evaluate the accuracy of our inner approximations to the 
open system vis-a-vis the closed system across a broad range of parameter values. In Section[5]we investigate the 
possibility of damped oscillatory behavior in the open system using two different approaches. We end with a brief 
discussion in Section|6] 



2 QSSA for the enzymatically open system 

In attempting to recapitulate standard QSSA analysis in the case of the open system ( [T0| ), we immediately face a 
significant obstacle. The usual first step is to exploit the conservation of total enzyme (free plus bound) to eliminate 
the differential equation for free enzyme concentration, but the inclusion of enzyme input and removal in the open 
system invalidates the balance relation ( [3b] ). Instead, the rate of change of total enzyme concentration equals the 
difference between the rate of enzyme input and the rate of enzyme removal: 



dE dC 

- + - =k -L. 3 E (11) 



This modified conservation relation is more usefully recast in integral terms: 



E(t) + C(t) =E + f (k 3 - k_ 3 E)dt (12) 
Jo 

The total concentration of enzyme in free and bound form is simply the initial enzyme concentration plus the 
sum over time of the net rate of enzyme accumulation. Though we are unable to reduce the dimensionality of the 
open system using this new balance relation, it does enable us to find closed form approximations to the enzyme 
and complex dynamics in the initial transient period. 
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2.1 Approximation to enzyme dynamics 

Solving (T2| for C and substituting into ( |10a| ), we obtain an integro- differential equation: 

d ^=k 3 -k- 3 E-k 1 ES+(k- 1 +k 2 ) (e + j\k 3 -k- 3 E)dt -E^j (13) 
Differentiating removes the integral and yields a second order differential equation 
d 2 E dE 

—Y + (k- 3 +kiS + k- X +k 2 )— + (k- 3 (k-i +k 2 ))E = k 3 (k-i + k 2 ) (14) 
dt at 

whose solution requires specification of initial conditions for E and ^ . 

Having S multiplying the first derivative term of equation |T4| renders it nonlinear and would appear to stymie 
attempts to find a simple closed form solution. At this juncture, we invoke the quasi-steady-state assumption: 
we reason that if S evolves on a much slower time scale than E, so that J| w 0, then we may treat S w So as a 
quasi-static parameter. (We investigate conditions under which this assumption is justified in Section|3]) 

Under this assumption, equation ( p4| becomes a linear second order ODE, which we proceed to solve by 
considering its homogeneous and inhomogeneous components in turn. The characteristic equation for the homo- 
geneous component 

r 2 + (£_ 3 +hS + k-i +k 2 )r+(k- 3 (k-x +k 2 )) = (15) 

has two roots: 



n, 2 = i ( -k- 3 -k x S -k-i-k 2 ± x / (£_ 3 + hS + k-i+k 2 ) 2 - 4k- 3 (k_i +k 2 ) 



\ k3-^-L,-^(L3^ + t, + Jl- ( k Jtt+Z 2 ^ 1 (,6) 



k-3+h{S + K M ) ( 4k- 3 K M 

V h{k-,/h+S + K M y 1 

The homogenous component of the solution is then Eh(t) — Ae nt +Be r2t , where A, B will be determined from 
boundary conditions. Note that both r\ and r 2 necessarily have negative real components (we consider the possi- 
bility that the roots r\ i2 are complex in a subsequent section). 

To find the inhomogeneous component of the solution, we make the ansatz Ei(t) — Ct 2 +Dt + F. Upon substi- 
tution into (14} , we find that Ei(t) = Ke, where Ke — k 3 jk- 3 is the enzyme accumulation constant for the reaction 
([9]). Note that Ke is not dimensionless, but is in fact in units of (enzyme) concentration. 

Thus, the full solution to (T4| ) is 

E{t) = E h {t)+Ei{t) —Ae nt +Be r2t +K E (18) 

For time t = 0, the initial conditions E(0) = Eo and ^ \ f=Q — k 3 — k- 3 Eo — k\EoSo determine A and B, 

-k 3 +E (k- 3 +kiSo)-r 2 (K E -Eo) 
A — (19a) 

r 2 -n 

B = k 3 - E (k- 3 -kiSp) + n (Ke - E ) 
r 2 -n 



which completes the inner solution for enzyme dynamics during the initial transient period. 
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2.2 Approximation to complex dynamics 

We find an approximate solution to the dynamics of the complex in a similar fashion. Substituting ( [T8} into ( |10b| ), 
we obtain 

^ + (k-x + k 2 )C = hS (K E +Ae r ^ +Be r2t ) (20) 



which is again linear if we treat S as constant. The homogenous solution to ( |20| ) is 

C h (t) = De-^- 1 ***)' = De~ klKut (21) 

where D is a coefficient to be solved for. Our ansatz for the inhomogeneous solution is Q(t) — Xe rit + Ye rit + Z, 
which we substitute into ( [20] ) and solve for X,Y,Z: 

x = — rrr* = /A r 5 ° (22a) 

7 = -7?V 50 = ^ ( 22b ) 

k 3 ki K E 
k-3(k-i+k 2 ) K M 

The full solution is then 

C(t) = C H (t)+Q(t) = So ( *f 1 ' + B6 \ + ^) +De- k ^ (23) 
V n /ki + ^ M r 2 /fci +Km K M J 

For the initial transient period, we use the t — initial condition C(0) = to fix D: 

D = _ 5o (_iL, + (24) 

\ri/ki+K M r 2 /ki+K M K M J 

The equation for C(f ) in the initial transient period becomes 

c ») - » « IT"' + B( ^; IT"' ) + ¥■ ('-'-**-') 

V n/ki+K M r 2 /ki+K M J K M V / 



Further simplification of p5| yields 



C(f) = 50 ( + ^TIm + Km) (26) 

The expressions we obtain for the approximate form of the enzyme and complex dynamics during the initial 
transient phase of the open reaction are more complicated than those for obtained by applying the QSSA for the 
closed system. Instead of a simple exponential relaxation to the quasi- steady- state, we find for the open system 
a double exponential modulation of the enzyme (respectively, complex) level set by the enzyme accumulation 
constant Ke. 

The substrate and product dynamics of the open system may be calculated by substituting fTS] ) and ( [26] ) into 
ODEs |T0b} and ( [TQd| ). 

^ = k-xC(t)-k x E(t)S (27a) 
at 

dP 

— = k 2 C{t) (27b) 



8 



Reznik, Segre, Sherwood 



For the remainder of the paper, we refer to the (non- autonomous differential-algebraic) system composed of 
Equations (T8) , ([26), ( |27a| , and p7b| as the approximate open system. We computationally investigate solutions 
to the approximate open system in comparison to solutions to the 'full' open system {TO) and the closed system 
in Section 0] 



2.3 Relation to classical Michaelis-Menten kinetics 



The closed system (|2) may be viewed as a special case of the open system ( [To) in which = k-3 = 0, which we 
may investigate using the above approach. In this case, Equation ( fl4) simplifies to 



d 2 E , ,dE d 2 E , ,dE 

+ (k x S + k-i+k 2 ) — = —2+k 1 (S + K M )- ] -=0 (28) 
dr dt dt z dt 

with initial conditions £(0) = Eq and ^ | Q = —kiEoSo. Again treating S as a parameter, the solution to (28 1 is 



E(t) =E - f l-e- k ^W) (29) 



Invoking the original conservation relation ([3b]), we immediately obtain an expression for the complex: 



C (t) = C 1 _ ,,-ti (s +K m ),\ m 

So + K m \ J 



Thus our approximations for the open system pass a basic consistency check: Equations ( [29] ) and ( [30) exactly 
match those obtained in classical derivations of Michaelis-Menten dynamics QOO. 



3 Time scales and quasi-steady states 

A key assumption underlying the derivations of the previous section is that the level of substrate does not decline 
appreciably from its initial value. This is obviously not true for all time, but as in other QSSA analyses, we expect 
the assumption to hold during an initial transient period as the levels of enzyme and complex evolve on a much 
faster time scale than the level of substrate. In the initial transient period, Equations (TS) and ( [26) should provide 
accurate inner approximations to the full dynamics of the reaction. The first aim of this section is to estimate the 
fast time scale of the open reaction fTO) , that is, to estimate the duration of the initial transient period. 

QSSA analyses of the closed reaction ([2) show that by the end of the transient period, the enzyme and complex 
reach a quasi-equilibrium with respect to one another. This quasi-equilibrium is maintained as the substrate level 
slowly declines and the reaction eventually extinguishes itself. In the open system, enzyme (and hence complex) 
levels are altered exogenously and at potentially rapid rates, so we do not necessarily expect the dynamics of the 
full system to be organized by slowly changing quasi-equilibria for enzyme and complex in that case. Nonetheless, 
we may derive expressions for the enzyme and complex quasi-equilibria that may be operative at the end of the 
transient period, and compare these expressions with those for the closed reaction. This is the second aim of this 
section. 

The section's third aim is to establish conditions on the reaction parameters under which we may be con- 
fident that our inner approximations accurately reflect the dynamics of the full open syste m fTO) . In particular, 
we seek inequalities in non-dimensionalized form which, when satisfied, imply equations (|18) and ( [26) closely 
approximate the true solutions E(t) and C(t) for the open system during the initial transient period. 
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3.1 Fast time scale (enzyme-complex) 

Initially the levels of enzyme and complex change quickly, evolving much more rapidly than the level of substrate. 
Thus we note first that the exponential form of equations fTS] ) and ( [26] ) present two immediate candidates for the 
fast time scale of the open reaction, namely the characteristic time scales r\ and r^. Since \ri\ > \n\, we expect 
the dominant fast time scale (and hence the time span of the initial transient period) to be 

= N = h(So + K M )+k- 3 (31) 
We observe that this expression differs only slightly from the usual estimate for the fast time scale in the closed 
system, which is 

tf= h(So+K M ) (32) 

In particular, the rate of enzyme removal k-3 reduces the duration of the transient period. The rate of enzyme 
input, by contrast, does not affect the estimate of tf. 

We have neglected the contribution of r\ in setting the time scale of the initial transient. Indeed, we argue that 
under most biologically relevant parameter regimes, fa\ ^> \r\\. Let 8 = (4k- ^Km)/ [^1(^-3/^1 +S + Km) 2 ], and 
reexpress the discriminant of \\5\ as y/i — 8 « 1 — 8/2 — &(8 2 ) for <5 1 (using the binomial rule). When this 
inequality holds, 

n ~ (33) 



n-n 



4 

2 



(2+2) ~-{k-3+h(S + K M )) (34) 

k- 3 +ki(S + K M ) { t t , $ S\ n . f „ 

2 \ l ~ l+ 2~2J = -\ k -3+ k n s + K M)) ~ n ( 35 ) 



Thus when 5 1, we have r\ w r2^/4, so that r\ <C r2, and we are justified in neglecting the contribution of r\ 
in setting the fast transient time scale. 

Two biologically plausible ways of satisfying the condition 8 1 are immediately apparent. First, when the 
level of substrate is relatively high (as compared to the product of the Michaelis-Menten constant and the enzyme 
removal rate), the denominator &-3/&1 +S + Km is large, and so 8 is small. Relatively high initial substrate levels 
are assumed in standard QSSA analyses, and over the course of the initial transient phase, we expect little change 
from the starting level of substrate (an assumption whose validity we investigate in Subsection |3 . 3 | below) . Second, 
when k-3 is small (relative to the Michaelis-Menten constant of the reaction), 8 will be small regardless of the 
substrate level. Enzyme degradation rates are typically at least two orders of magnitude lower than reaction rates 
for complex and product formation; thus for biologically relevant parameter regimes we may reasonably assume 
k-3 Km and hence 8 <C 1. 



3.2 Quasi-steady states 

In the closed system, the relative levels of enzyme and complex quickly equilibrate with respect to the level 
of substrate. This rapid equilibration occurs within the initial transient period, and subsequently the enzyme and 
complex quasi-equilibria evolve on a slow time scale that tracks the slow depletion of substrate. Though we do not, 
in general, expect the enzyme and complex dynamics in the open system to be as strictly governed by quasi- steady 
states after the initial transient period, we may calculate the quasi-steady states just as for the closed reaction. 



We begin by observing that the dominant characteristic decay time for equations ( [18] ) and (26) is tf. At the 
end of the initial transient period, when €{$f) time has elapsed, those terms multiplied by £ r2 cipproach zero, 
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while (assuming 8 <C 1 as above) e nt remains near unity. Hence for E, the quasi-steady state of the enzyme 
concentration, we find that at the end of the initial transient period 



A + K E 

-k 3 + k- 3 E + k\E S - 



r 2 (k 3 /k- 3 -E ) k 3 



-k 3 +k- 3 Eo + kiEoS 

-(fc_i+fc 2 + fciS + fc-3) 
EqKm + k 3 /k\ 
S + K M + k- 3 /ki 



+ ■ 



(36) 
(37) 
(38) 



The above expression bears close resemblance to the quasi-steady state one computes for the closed reaction, 
E — EoKm/(So + Km), with two intuitive differences. First, the putative quasi-steady state for the open reaction is 
elevated by k 3 jk\ , the ratio of the enzyme influx to the enzyme-substrate binding rate, a quantity representing the 
net accumulation of free enzyme at equilibrium. Second, the putative quasi-steady state is reduced by k- 3 jk\, the 
ratio of enzyme removal to the enzyme- substrate binding rate, which reflects the net subtraction of free enzyme at 
equilibrium. Whether the putative quasi-steady state for the open system is higher or lower than that of the closed 
system depends on the relative magnitudes of k 3 and k- 3 . 

Our derivation of C, the quasi-equilibrium for complex concentration at the end of the initial transient period, 
proceeds similarly: 



CKk lS ( , A , + ^-) (39) 

\r\ +K-\ + k 2 k- \ + k 2 / 

= k s ( -k 3 +k- 3 Eo + kiEpS E -k 3 /k- 3 k 3 /k- 3 \ 

~ l -Ik i+k 2 +hS)(k x +k 2 ) (*-i+*2+*i*)(*-i+*2) + k x + kn _ ^ ) 

\ k_i+k 2 +kiS+k_ 3 / 

_ 7 v f(Eo-k 3 /k-3)(k-i +k 2 )-kik 3 S/k- 3 + k 3 /k- 3 \ (41) 



(k-i +k 2 +k x S)(k-i +k 2 ) k-i +k 2 

Eq \ EqS 



k-i -\-k 2 -\-k\S J S + Km 



(42) 



Rather remarkably, this final expression is identical to the quasi-steady state for the closed reaction ([5]). 

Our computation of the quasi-steady states relies on two assumptions which may not necessarily hold true for 
the open reaction, depending upon parameter choices. First, we explicitly assumed that free enzyme and complex 
levels quickly equilibrate relative to one another. If free enzyme is added to or removed from the open system 
at a rate comparable to the usual enzyme-complex equilibration rate, our assumption does not necessarily hold, 
in which case the quasi-steady states calculated above are not meaningful for any substantial period. Second, we 
implicitly assumed that the amount of total enzyme (free and bound) remains essentially constant in the post- 
transient period. For relatively small values of k 3 and k- 3 , this assumption may hold nearly true for an extended 
period, during which time expressions ( [38] ) and ( [42} may provide reasonably accurate approximations to the 
enzyme and complex concentrations of the system. However, even for relatively small values of k 3 and k- 3 , the 
accumulating deviation from constant total enzyme level will not remain negligible indefinitely. Thus ( [38] ) and 
([42} may not accurately estimate enzyme and complex levels in the long term. Over shorter periods, however, and 
in the transient period in particular, we may be assured that the approximations from this and prior sections are 
valid, providing certain conditions are met, as discussed below. 
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3.3 Validity conditions 



To ensure the accuracy of the inner solutions (T5] ) and ( [26] ), we require foremost that the substrate concentration 
not deviate significantly from its initial level during the initial transient period. That is, for time t < @(tf), we 
require S(t) w So, or equivalently, |AS/So| <C 1. 

During the initial transient period, the total change in substrate concentration AS is coarsely bounded by 
\Smax\ • tf, the product of the maximum rate of change in substrate and the length of the transient period. Thus 
using the formula 



AS 


1 


dS 






W 5b" 


dt 


max 



(43) 



we certainly avoid underestimating the relative change in substrate. In the QSSA for the closed system, it can be 
shown that S = k-\C(S,t) — k\E(S,t)S achieves its maximum at (S,t) ~ (Sb,0). F° r small values of k 3 and k- 3 , 
this is nearly true in the open system as well (our numerical investigations in Section [4] demonstrate that this is not 
always the case, e.g. when k-3 is small and k 3 is relatively large), and using this approximation we then find 



AS 



kiEoSo 



Eo 



S (k-i+k 2 +kiS + k- 3 ) 
At this point we define several dimensionless parameters: 



S+K M + k-3/h 



K m 



k-x 
k 2 



k\K m 



(44) 



(45) 



Here a and p augment the parameter set (a, 77, k) familiar from the canonical nondimensionalization of the 
closed system reaction [ 5 ]. These dimensionless parameters may be used in conjunction with the rescaled variables 
e — E/E, c = C/C, s = S/So, 1 — t /tf to nondimensionalize the equations for the open system ( 10). 
Rewriting |44|) in dimensionless form, the condition |^5/5q| <C 1 becomes 



7] < 1 + a + a 



(46) 



When this inequality holds, we may be confident of the accuracy of our inner solutions during the initial transient 
period. We note that ( [46] ) is weaker than the corresponding inequality for the closed system, namely 77 <^ 1 + o 
0. (Furthermore, the two inequalities coincide as k-3 goes to zero, that is, for negligible enzyme removal.) 
Thus in any parameter regime for which the standard QSSA is valid, equations (T8| and ( [26] ) hold (with similar 
approximation errors) for t < @(tf). Numerical simulations support this conclusion, and show additionally that 
within a broad range of parameter values the approximations remain accurate well beyond time tf. 



4 Numerical results 

Subject to some constraints on rate constants, we expect the approximations derived in Sections [ZT] and |Z2] to 
provide good estimates of the enzyme and complex levels of the open system during the transient phase. It is not 
entirely clear from the derivation just how long the key assumption, namely that S w So, remains valid, nor how 
the domain of validity is affected by the rates of enzyme addition and removal, especially when one or both of 
&3,&_3 approach the magnitude of the other reaction parameters. To gauge the efficacy of our approximations for 
different rates of enzyme input and removal, we numerically computed solutions for the full set of ODEs fTQ| ) of 
the open system and for the approximate open system |T8| , ( [26] ), p7| for broad ranges of values for the rates of 
enzyme addition and removal (parameters k 3 ,k-3, respectively). For the approximate open system, a differential- 
algebraic system of equations, we used the analytic expressions derived for E(t) and C(t) in the transient phase 
(equations (T8| and ( [26] ), respectively), and we numerically solved the (nonautonomous) differential equations 
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|27a| and ([27b]), in which E(t) is given by equation ( 18 ) and C(t) is given by equation \26) . We also computed 
solutions to the closed system ([2} in order to better assess the utility of our approximation. Numerical integration 
was done in MATLAB using the ode 15s stiff integrator with relative and absolute error tolerances of 10 -10 . We 
consider two regimes of initial conditions, the standard QSSA regime with So ^> Eo, and the reverse QSSA regime 
with Eo^> So. The simulated initial conditions were (So, Eo) — ( 1 , 0. 1 ) and (So, Eo) = (0. 1 , 1 ) for the standard and 
reverse QSSA regimes, respectively. 

In each regime of initial conditions, we perform a set of one-dimensional parameter sweeps: For each value 
K e {10 -4 ,0.1, 1} (a sampling of parameter values which is representative of the range of different scenarios of 
enzyme addition and removal), we fix a parameter k^i G {—3,3} at ki = K, and we simulate the open system 
(both the full set of ODEs and our approximation) for 1600 values of k-f logarithmically spaced in the interval 
[0.0001 , 5] . We display the portion of each simulation run spanning the time interval [0, 1] (tf) <C 1 for all parameter 
combinations considered), and we note that tf (the fast complex time scale) is less than 1 for every pair of £3, k-3 
values used. For the results presented below, we fixed the other parameters (which determine Km) at the same 
standard values in each simulation, namely k\ — k-\ — Jc2 = 1. (We also performed the same computations using 
other values for &i , &_i , &2 in order to sample a range of Km values, and we obtained results very similar to those 
presented below.) We did not exhaustively search the parameter space, but the selected examples we show give 
a comprehensive picture of the behavior of the full open system and our approximation under a wide range of 
conditions, including biologically relevant parameter ranges. 



Solution manifolds We present our numerical results in Figures[T[fl"2|below, with each figure displaying the results 
of one parameter sweep (1600 individual simulations). We first organize the figures according to initial conditions 
(standard followed by reverse QSSA regimes), and then for each regime of initial conditions, we group the figures 
according to which of £3, £-3 is the 'free' parameter being varied in the sweep (first varying k-3, then varying £3). 
A solution to either initial value problem {To} or fT8} , ( [26] ), ( [27} for a fixed value of the free parameter traces a one- 
dimensional curve in the four-dimensional phase space (E,C, S, P), and if we smoothly vary the free parameter, we 
obtain a two-dimensional manifold of solutions in the five-dimensional space (E,C, S, P, ki), where ki E {^3,^-3} is 
the free parameter. The results of a parameter sweep provide an approximation to these manifolds, and each figure 
depicts the solution manifolds for (To| and ( [T8} , ( [26] ), ( [27} for a single parameter sweep. Each figure comprises 
five subplots, and each of the first three subplots shows the projections of the solution manifolds onto one of the 
subspaces (E,S,ki), (C,S,ki), and (E,C,ki). These projections may be thought of as the union, taken over the set 
of values swept for the free parameter k, of the two-dimensional phase plane diagrams of the solutions obtained at 
fixed values of ki, for the (E,S), (C,S), and (E,C) planes, respectively. We omit projections involving the phase 
variable P, since the equations for product formation decouple in ([2}, (T0| and (T8| , ( [26] ), ( [27} , and thus this 
variable provides essentially redundant information. 

Each three-dimensional subplot depicts three surfaces: the solution manifold of the full open system ( [To} , 
denoted «y#o, in red; the solution manifold for the approximate open system fTo} and ([18] [26] [27}, denoted 
in blue; and the manifold of the full solution to the closed system ([2), denoted «^c> m greerjj 

We may think of the closed system ([2} as an approximation to the open system (To} that presumes the effects 
of enzyme input and removal to be negligible. For example, if £3, £-3 <C £1,^2? k-i, input and removal occur on 
a time-scale much slower than the individual chemical reactions which constitute the conversion of substrate to 
product. Viewed in this light, the projections of the manifold m Figures [T||T2] illustrate the degree to which 
this assumption of negligibility may mislead. Note that there are only two versions of «^c, one for the standard 
QSSA regime and one for the reverse QSSA regime, so that the <Jtc manifolds shown in Figures [IJj6] are identical 
to one another, as are the manifolds shown in Figures |7}{T2| 

1 Neither A3 nor k-3 appear in {5}, so the solution manifold is simply [0.01, 5] x where ^ denotes the solution of the initial 
value problem |2] over the time interval [0, 1] (with &i = = &2 = 1, and with initial conditions corresponding to the either the 
standard or reverse QSSA regimes) projected onto the appropriate phase plane. 
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In each three-dimensional subplot of Figures [T}fT2] the value of the free parameter is measured on a log- 
arithmic scale on the vertical axis, and two specific phase space variables are measured on the two horizon- 
tal axes. Time is not measured explicitly as a coordinate in these phase plots, but the temporal orientations of 
the solution manifolds may be inferred by noting the location of the manifolds' edges corresponding to the 
initial conditions. For the standard QSSA regime, the solution manifolds "start" at the edge corresponding to 
(E,C,S,P,ki) — (0.1,0, 1,0,&;), and for the reverse QSSA regime, the solution manifolds start at the edge corre- 
sponding to (E,C,S,P,ki) — (1,0,0.1,0,/:/). The end of the transient period, which occurs at time t « tf < 1, is 
indicated by a line drawn on each of the solution manifolds. The standard estimate for tf in the closed system, 
given by ( [32} , is drawn as a dashed line, and the estimate for tf derived in this paper for the open system, given by 
( [3T] ), is drawn as a solid line. For visual ease, both the dashed and solid lines are drawn on the <Jto (red) and 
(blue) manifolds in black, and on the (green) manifold in magenta. 

Approximation error We have mentioned repeatedly that QSSA-based analyses presume S ^ So during the pre- 
steady-state transient. The fractional decrease in substrate during this period, 8s = (S(tf) — So) /So, provides one 
measure of the error for the QSSA approximation. As discussed in (5), given a prescribed tolerance for 8s, one 
useful heuristic for gauging whether the QSSA can be used (for the closed system) is the condition 

E ° ~ 71 <S S (47) 



K M + So 1 + ct 
For the open reaction, this condition becomes 

E _ V 



< 8 S (48) 
K M + S + k- 3 /ki l + a + a 

Condition ( [48] ) is weaker for the open system than condition ( [47] ) is for the closed system; both conditions imply in 
particular that the approximations should be accurate in the standard QSSA regime. Furthermore, we may expect 
the inner approximation for the open system to be more accurate at higher rates of enzyme removal. 

For each parameter sweep we present, we also plot a different measure of the approximation error. The fourth 
(two-dimensional) subplot in each of Figures [Tf|T2] records the signed difference between the actual solution to 
(T0| and the approximate solution fT8]|26]|27} . The fifth (two-dimensional) subplot in each of Figures[T[]T2]records 
the signed difference between the solution to the full open system ( [To| and the solution of the full closed system 
{2}. In both cases, the signed error is taken at the end of the transient period, time tf, where this time is given by 
the new estimate ( [3T] ): 

E x (tf)=X(tf)-X(tf) (49) 



Here X is one of E, C, S; X(tf) denotes the X component of the actual solution to ( 10 ) at time tf, mdX(tf) denotes 



the X component either of the approximation ( 18 ), ( [26} , ( 27 ) or the full closed system at time tf. The signed error is 
calculated for each simulation as the free parameter kt G {£3, k-3} is varied. We use the signed difference between 
the actual and approximate solutions, rather than the absolute value of their difference, in order to make clear 
when the approximation underestimates or overestimates the true values of E,C,S. 



4.1 Standard QSSA regime 

When So^> Eo, our validity condition for inner approximation {T8}, j26} is satisfied by a substantial margin. As 
reflected in the proximity of the red and blue manifolds in Figures [Ifj6] for a broad range of parameters our inner 
approximation provides a highly accurate estimate of the behavior of the open system up to, and often well past, 
times @(tf). Furthermore, the closed system poorly approximates the behavior of the open system, especially for 
non-negligible levels of enzyme input. 
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Varying enzyme removal rate The limitations of the closed system for approximating the open system are quite 
apparent in simulations where k-3 varies, shown in Figures [T}j3] When the rates of enzyme input and enzyme 
removal are both very low, the closed system does match the behavior of the open system well: the manifolds .Jto 
(red) and Jtc (green), as well as ^a (blue), all lie very close together for values of k-3 two orders of magnitude 
and more below the value of Km, as seen in Figure [T] Above that rate of enzyme removal, the behavior of the open 
and closed systems begin to deviate substantially, with the open system overestimating the level of free enzyme, as 
is to be expected. Total enzyme levels decline as free enzyme is removed from the system relatively rapidly, rather 
than remaining approximately constant. Specifically, the results recorded in Figure [Te] indicate that the signed 
error between the open and closed systems at time tf is <^(10 -4 ) for k-3 — 10 -3 and increases to <^(10 -2 ) as k-3 
approaches 1, with the sign of the error in the E and C components always positive and that of the S component 
always negative. 

At more rapid rates of enzyme input, the shape of <Jto differs greatly from that of «y#c, indicating that the 
closed system poorly approximates the behavior of the closed system at almost any level of enzyme removal. 
This can be seen in Figure [2] in which the the rate of enzyme input is an order of magnitude lower than the other 
reaction rate constants. Here <J£o and intersect for at a unique value of &_3(from Figure |2e] at k-3 « 1); at 
this parameter combination, the total amount of bound and unbound enzyme in the open system remains roughly 
constant over the transient period, as in the closed system. Elsewhere, the ratio between £3 and k-3 determines 
whether the closed system over- or under-estimates reactant levels for the open reaction. The position of 
vis-a-vis Jto reflects underestimation of enzyme and complex levels and overestimation of substrate levels when 
enzyme input is high and enzyme removal is low. Enzyme accumulates in free and bound form, rather than re- 
maining constant. When enzyme input is low and enzyme removal is high, the manifold configuration is reversed. 
In both cases, the numerical results confirm our intuition regarding the behavior of the closed system considered 
as an approximation to the open reaction. In Figure[3] which depicts results when the rate of enzyme input matches 
the other reaction rates, we see that for all values of k-3, and .Jtc intersect only at (Eq,Cq,Sq). The closed 
system always underestimates enzyme and complex levels and overestimates substrate levels (Figure [3e|. The 
range of absolute discrepancies between the full open and closed systems grows from <^(10 -3 ) to <^(10 - ~as k-3 
increases from 0.0001 to 1 (Figures [Te}|3e] ). 

In each of Figures [TJj3] the manifolds ^a and <Jto lie very close to one another for all values of k-3, indicating 
that the inner approximation to the open system is very accurate, even when the rate of enzyme input is very rapid. 
In particular, the magnitude of the signed error between the full system and the inner approximation at tf recorded 



in Figures Id - 3d is small, ranging only from <^(10 ) when k-3 — 0.0001 to <^(10 ) when k-3 — 1. The position 
of the manifold relative to the manifold reflects consistent underestimation of the level of enzyme and 
overestimation of the levels of complex and substrate. Since the inner approximation was derived by assuming 
constant substrate, the magnitude of r^, the dominant, negative exponent for enzyme and complex dynamics in 
Equations fTS] ) and ( [23] }, does not decline as the substrate level falls. The decaying exponential component of 
Equation (|18|) is overestimated, so the predicted enzyme level is lower than the actual enzyme concentration 
obtained from the full set of ODEs ( [T0| ) of the open system. Furthermore, the approximate complex level given 
in Equati on p3| is directly proportional to So, and hence is overestimated. The substrate level, obtained via 
Equation ( |27a| ), declines at a rate proportional to the enzyme level, and so it is underestimated in tandem with 
the overestimation of enzyme concentration. We note that the magnitudes of the signed errors in E,S,C decrease 
as k-3 grows, particularly the signed errors in enzyme and complex levels. This accords with the behavior of the 
error condition ( [48] ) as k-3 increases, discussed above. 

The curvature of the «y#o and <J^a manifolds seen in each of Figures [T}|3] at higher levels of k-3 also accords 
with our intuition. For a given rate of enzyme addition, £3, we expect higher rates of enzyme removal to lead to 
lower levels of free enzyme; hence the E-C and E-S projections of <Jto and J£a fold in the direction of lower en- 
zyme levels as k-3 increases. Furthermore, more rapid enzyme depletion reduces the amount of enzyme available 
to catalyze the conversion of substrate to product, and thus at larger values of k-3 the level of substrate remaining 
at tf is greater. 
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(a) Substrate-Enzyme (b) Complex-Enzyme (c) Complex-Substrate 




(d) Signed error at tf, full open (e) Signed error at tf, full open 
system vs. open approximation system vs. full closed system 

Fig. 1: Phase plane portraits and signed errors in the standard QSSA regime as k-3 varies, with £3 = 0.0001. 



Varying enzyme input rate Figures [4j|6] show the solution manifolds «y#o, and «y#c when £3 is varied and 
k-3 is held fixed. Qualitatively, the respective solution manifolds at different k-3 values differ only little from one 
another; the chief difference between them is the change in scale. As in the case above, where £3 is held fixed 
and k-3 varied, we find that the open approximation outperforms the closed system as an estimator of the true 
behavior of the open reaction. The degree to which the open approximation is superior increases with higher k-3 
values, as one would expect. 

For each fixed value of k-3 the <Jto and manifolds deviate substantially from the manifold once 
&3 is above a certain value, which depends on k-3. In each of Figures |4j|6] the Jto and manifolds form 
large 'flaps' that extend far away from the <y#c manifold in the directions of increased enzyme and complex and 
decreased substrate. These flaps correspond to parameter regimes in which the rate of enzyme input is relatively 
rapid and thus assuming that total enzyme remains constant is a clearly poor approximation to the behavior of 
the open reaction. Free enzyme accumulates, as does complex, and substrate is rapidly degraded. The deviation 
between the <Jto~^A manifolds and the Jtc manifold is obviously greatest at the uppermost rates of enzyme input, 
where £3 significantly exceeds k-3 and begins to attain the same order of magnitude as Km- We note, however, 
that even at these high values of £3, the open approximation (T5] ), ( [26] ) lies remarkably close to the solution of the 
full equations for the open system up to time tf. As the value of £3 decreases, the discrepancy between the open 
and closed systems decreases, as expected. As shown in Figures |4e}{6e| while for each k-3 value examined the 
signed error between the full open and closed systems rises to <^(1) as £3 increases beyond approximately 10 _1 , 
within the three lowest decades of £3 values, the discrepancy between the full open and closed systems ranges 
from <^(10 -4 ) when k-3 = 0.0001 to <^(10 -2 ) when k-3 — 1. As an estimate of the full open reaction, the open 
approximation is always at least as accurate as the closed system, even at low levels of enzyme input, and it is 
typically more accurate by more than an order of magnitude. For low £3 values, this is difficult to discern in the 
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(a) Substrate-Enzyme (b) Complex-Enzyme (c) Complex-Substrate 




(d) Signed error at tf, full open (e) Signed error at tf, full open 
system vs. open approximation system vs. full closed system 



Fig. 2: Phase plane portraits and signed error in the standard QSSA regime as k-3 varies, with £3 =0.1. 



three-dimensional subplots of Figures |5J|6] due to the size of the flaps of the and manifolds. Figures pd]- 
[6d| however, show that for each £3 value the signed error between the full open system and the open approximation 
is <^(10 -4 ) over the three lowest decades of k-3 values and rises to at most <^(10 -2 ) for k-3 =0.1. Comparing 
with Figures |4e||6e| we see that at identical parameter values, the open approximation is generally one to two 
orders of magnitude more accurate than the closed system as an estimate of the full open system. 

The position of the <J^a manifold relative to that of the .Jto manifold, at all £3 values, indicates that the 
open approximation consistently underestimates the level of free enzyme and overestimates the level of complex 
in the open reaction. This situation can also be seen in the subplots of the signed approximation error, [4d|j6d] 
which also show that the estimate for substrate levels is quite accurate at nearly all £3 values. The reasons for 
this arrangement of estimation errors are largely the same as those given above. As an addendum to the previous 
explanation, we note that the open approximation depends in part on the rate of addition of enzyme to the system 
being small relative to other reaction rates. In particular, we assume in Equation ( [44] ) that the maximal rate of 
change in substrate concentration occurs at time t — 0, which is certainly the case in the closed reaction and is 
valid at lower values of £3. As £3 increases, this assumption becomes less accurate; to consider a limiting case, 
the introduction of enzyme via an instantaneous pulse will produce an immediate, substantial increase in the rate 
of substrate degradation, and thus a greater deviation from the presumed value S ^ So which appears in (T5] ), p6| ). 
We note also that the error condition ( [48} remains unchanged as £3 varies. Hence we would expect the increase in 
signed error with increasing £3 shown in the Figures |4d|{6d| which is the reverse of what is seen when £3 is held 
constant and k-3 is varied. 

Before turning to the results for the reverse QSSA regime, we make two related observations which may 
aid in interpreting the figures in this subsection and the next. First we note that the intersection of the and 
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(a) Substrate-Enzyme (b) Complex-Enzyme (c) Complex-Substrate 



(d) Signed error at tf, full open (e) Signed error at tf, full open 
system vs. open approximation system vs. full closed system 

Fig. 3: Phase plane portraits and signed error in the standard QSSA regime as k-3 varies, with £3 = 1. 



«y#c manifolds marks the location of parameter combinations for which the total amount of enzyme remains 
approximately constant in the open system over the transient period. This occurs roughly where £3 « k-3, as is 
somewhat more easily seen in Figures [TJj3] than in Figures [4}j6] To reiterate the conclusions detailed above, when 
&3 < k-3, the closed system overestimates the level of enzyme, and when £3 > k-3, the situation is reversed. 

Second, loosely speaking, each of Figures [TJj3] is "identical" up to an extension (scaling) and rotation of the 
<y#c and manifolds. To see this relationship, observe first in Figure [Ta| how the <y## and <J£a manifolds lie 
close to the <Jtc manifold at low values of k-3 and peel away from it as k-3 increases, folding over in the direction 
of lower enzyme concentration. The same folding is seen in Figure |2a| but here the fold intersects Jtc where 
k3 « k-3. The <y## and J%a manifolds extend substantially past the «y#c manifold in the directions of higher 
enzyme concentration and lower substrate concentration; this is most noticeable in the lower portions of <Jto and 
^A- This extension relative to the unchanging <J£q manifold is due to the input of free enzyme, which increases 
enzyme concentration and accelerates depletion of substrate. The increase in £3 between Figure [la] and Figure |2a| 
effectively rotates and stretches <Jto and ^a relative to Mc- At the highest level of £3, shown in Figure [3a] <y#o 
and ^a retain their fold shape, but their elongation is much more pronounced and the rotation is even greater, so 
that and ^a lie entirely on the other side of «y#c- 

Similar relations hold between the manifolds depicted in Figures [Tb] - [3b] and Figures [TcU3c| F urthermore, 
Figures [TJj3]in this subsection are related to one another in an analogous fashion, as are Figures |10H12| and Figures 
[7}|9]in the following subsection. 
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(a) Substrate-Enzyme 



(b) Complex-Enzyme 



(c) Complex-Substrate 




(d) Signed error at tf, full open (e) Signed error at tf, full open 
system vs. open approximation system vs. full closed system 

Fig. 4: Phase plane portraits and signed error in the standard QSSA regime as £3 varies, with k-3 — 0.0001. 



4.2 Reverse QSSA regime 

When Eo^> So, the validity condition |46| does not hold, i.e. 77 /(l + a + a) 1, and we cannot be confident that 
the open approximation p8| ), ( [26] ) will accurately capture the behavior of the full open system. Also uncertain 
is the accuracy of the closed system as an approximation to the full open reaction. The results of numerical 
simulations in the reverse QSSA regime, however, bear striking overall resemblance to the corresponding results 
for the standard QSSA regime. In particular, our numerical results indicate that the open approximation provides 
a quite accurate approximation to the full open system. The closed system approximates the open reaction well 
only when the ratio of £3 and k-3 implies that the change in enzyme level up to time tf is miniscule. 

The reverse QSSA regime results shown in Figures[7|j9] in which k-3 varies, and in Figures [l0j|T2] in which £3 
varies, match up very closely, at least qualitatively, with the corresponding results for the standard QSSA regime. 
As in the standard QSSA regime, and for the same reasons, the open approximation consistently underestimates 
the concentration of free enzyme and overestimates the concentration of complex and substrate, at least up to time 
tf, the estimated end of the transient period. The closed system over- or under-estimates enzyme, complex, and 
substrate levels, depending on the ratio of £3 and k-3. As can be seen in Figures [7df|12d| the difference between 
the true quantities in the full open system and the estimates from p8| , ([26) follows the same pattern in the reverse 
QSSA regime as in the standard QSSA regime when £3 and k-3 are varied. We point out, however, that past 
time tf the open approximation may begin to underestimate the level of substrate. This phenomenon can be seen 
in Figures [7c}j9c] where the C-S projection of the <J^a manifold intersects and folds behind the C-S projection 
of the manifold. This occurs in scenarios where k-3 significantly exceeds £3, such that the amount of total 
enzyme available, and hence the effective substrate degradation capacity, drops quickly and consequentially over 
the course of the reaction's early stages. We also note that in the reverse QSSA regime the closed system's estimate 
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(a) Substrate-Enzyme 



(b) Complex-Enzyme 



(c) Complex-Substrate 




(d) Signed error atZy, full open (e) Signed error at tf, full open 
system vs. open approximation system vs. full closed system 

Fig. 5: Phase plane portraits and signed error in the standard QSSA. regime as varies, with k—3 — 0.1. 



of complex and substrate levels roughly matches its accuracy in the standard QSSA regime, as can be seen by 
comparing Figures [7ef| 1 2e| with Figures [Te}|6e| 

The chief difference between corresponding sets of results from the standard and reverse QSSA regimes is 
quantitative. The signed error between the full open system and the open approximation in each of enzyme, 
complex, and susbstrate, measured at time tf, is typically four to ten times greater in the reverse QSSA regime 
than in the standard QSSA regime, at least when the free parameter is smaller than Km by an order of magnitude 
or more. To explain this phenomenon, we refer to condition ( |46| . In the reverse QSSA regime, 7] > cr, so for the 
condition to hold, a, the non-dimensionalized surrogate for k-3, must be quite large. At higher k-3 values, the 
accuracy of the open approximation should improve, and this is precisely what may be seen by comparing, for 
example, the error scale in Figure [T0d| (/:_ 3 = 0.0001) with that of Figure [T2d| (/:_ 3 = 1), or by examining any of 
Figures |7d}|9d| In general, we would expect the magnitude of error between the full open system and the open 
approximation at identical values of &_3and £3 to be higher in the reverse QSSA regime than in the standard QSSA 
regime, which is indeed the outcome we find numerically. 
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(a) Substrate-Enzyme 



(b) Complex-Enzyme 



(c) Complex-Substrate 




(d) Signed error atZy, full open (e) Signed error at tf, full open 
system vs. open approximation system vs. full closed system 

Fig. 6: Phase plane portraits and signed error in the standard QSSA regime as £3 varies, with k- 




Fig. 7: Phase plane portraits and signed error in the reverse QSSA regime as k-3 varies, with £3 = 0.0001. 
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(a) Substrate-Enzyme (b) Complex-Enzyme (c) Complex-Substrate 



I 

1° 




(d) Signed error at tf, full open (e) Signed error at tf, full open 
system vs. open approximation system vs. full closed system 

Fig. 8: Phase plane portraits and signed error in the reverse QSSA regime as k-3 varies, with £3 =0.1. 




Fig. 9: Phase plane portraits and signed error in the reverse QSSA regime as k-3 varies, with £3 = 1. 
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(a) Substrate-Enzyme 



(b) Complex-Enzyme 



(c) Complex-Substrate 



3 



(d) Signed error at tf, full open (e) Signed error at tf, full open 
system vs. open approximation system vs. full closed system 



Fig. 10: Phase plane portraits and signed error in the reverse QSSA regime as varies, with k-3 = 0.0001. 




(a) Substrate-Enzyme (b) Complex-Enzyme (c) Complex-Substrate 




(d) Signed error at tf, full open (e) Signed error at tf, full open 
system vs. open approximation system vs. full closed system 



Fig. 11: Phase plane portraits and signed error in the reverse QSSA regime as £3 varies, with k-3 = 0.1. 
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Fig. 12: Phase plane portraits and signed error in the reverse QSSA regime as £3 varies, with k-3 = 1. 
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5 Decaying oscillations? 

Neither the closed system ^ nor the open system (TU| ) is capable of sustained oscillations without substrate input. 
The reaction proceeds until all substrate is irreversibly converted into product. It is conceivable, however, that 
the approach to the steady state (£,C,S,P) — (£b 5 0,0,So) = E* could involve decaying oscillations, i.e. the fixed 
point F* is a stable spiral. 

It is a simple matter to linearize ([2) and (l0| ) about F* and calculate the four corresponding eigenvalues 
Ai , . . . , A4. For the closed system ([2), Ai = A2 = 0, and for the open system (To| , Ai = k- 3 , A2 = 0. The remaining 
eigenvalues for both systems are given by 



^3/ 



-(E ki +k-i+k 2 )± yj {Eoki +k-i +k 2 ) 2 -4E k l k 2 



(50) 



For decaying oscillations we require a negative discriminant, i.e. 



2E h k 2 > Eq k\ + k 2 _ l + k 2 + 2E ki k- 1 + 2k- 1 k 2 (5 1 a) 

k] k 2 1 k 2 k-\ k-\ 
k 2 E k\k 2 E k\ k 2 E k 2 



Since each of the terms in the right hand side of ( |51b| ) is nonnegative, satisfying that inequality entails satisfying 

2>E^ + p- (52) 
k 2 E k\ 



Writing A = Eo ^ , we see that satisfying inequality ( 52 > is equivalent to finding positive A such that > A 1 — 2A + 
1, which is impossible. Hence F* cannot be a stable spiral, and the system does not exhibit decaying oscillations. 

It has been shown elsewhere that for some systems, the reduced equations obtained via QSSA analyses may 
produce dynamical behavior at odds with that of the original system C~T][T8). ^ n Particular, the reduced system 
may be incapable of oscillation though the original system does in fact oscillate fTTl . Here we briefly investi- 
gate whether the QSSA analysis presented above may mislead in the opposite direction, that is, indicating the 
possibility of (decaying) oscillations when they are in fact impossible in the full system. 

Our point of departure is the observation that our method of solution for |T4| ) presented in Section ([2|, namely 
rewriting the first order nonlinear ODE ( |10a| ) as a (pseudo-) linear second order ODE |T4| , we obtain the equational 
form for a forced damped oscillator: 

^2 + 2y- + co 2 x = F(t) (53) 

In the case of a mechanical oscillator, such as a mass-spring system, 7 > represents the coefficient of friction, 
(0 represents the natural frequency of oscillation, and F(t) is the forcing term. In the enzyme dynamics equation 
(T4| for the open system, these coefficients are 

7= \{k^+kiSo + k-i+k 2 ) = l -{k^+ki{S + K M )) (54a) 

co 2 = k- 3 (k-i +k 2 ) = k- 3 kiK M (54b) 
F(t) = k 3 (ki+k 2 ) = k 3 kiK M (54c) 

Forcing is constant and proportional to the rate of enzyme input, the putative oscillation frequency is proportional 
to the rate of enzyme removal, and damping is (in the standard QSSA regime) dominated by the initial substrate 
concentration or Km- Viewed as a forced damped oscillator, the system would appear to have an natural frequency 
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which is independent of the initial enzyme concentration, a conclusion already at odds with the calculation above. 
(Obviously, any oscillatory solution may only be valid within the time frame for which Equation ( p4| provides a 
reasonably accurate approximation to the dynamics of the full system, that is, for time t < tf. If tf is less than the 
expected period of oscillation, then we may not be fully confident of realizing damped oscillations in the system.) 

Considering the approximate solution p8| we obtain for the second order ODE for enzyme, we would antici- 
pate the possibility of damped oscillations if r\ , r2 are complex valued, that is if 

y 4k- 3 K M 

*i(W*i+* + *i#) 2 P j 

or, equivalently, 



2y/KMk-3/k! - (K M + k- 3 /h) > So (55b) 



For any combination of positive parameter values, the maximum value of the lefthand side of inequality ( |55b| ) is 
zero, so that the inequality cannot be satisfied for positive initial substrate. Thus r\ , r2 are always real valued, and 
the QSSA approximation indicates that the system should not oscillate for physically admissible parameter values. 

We may also argue against the possibility of oscillations another way. Turning back to equations ([53])— ([54]) , 
the most propitious situation for oscillation is when there is no forcing, i.e. k 3 — 0, no enzyme input. In this case, 
decaying oscillations may potentially occur only if the system is underdamped, that is, when co > y. This requires 

1 



y/k- 3 kiK M > -(k-3+hiSo + KM)) (56a) 

4k- 3 k { K M > k- 3 2 + 2^3^ + K M )+k 2 (S + K M ) 2 (56b) 
> (k-3+hSo) 2 + 2k l {k l S -k- 3 )K M + kiK M 2 (56c) 



Condition (56c) implies that Km must lie between the roots of the quadratic, ^ [—(k\So — k-3) ±2y/kiSok-$\ . 
The requirement that Km be positive necessitates 



2^/k\Sok-3 > k\So - k-3 (57a) 

2> t^ + ft (57b) 

k-3 k\So 



Applying the same reasoning as we did for inequality ( [52] ), we see that inequality ( |57b| cannot be satisfied, and 
damped oscillations are impossible. 

In the vernacular of harmonic oscillators, thus we find that the enzymatically open system is always over- 
damped. Decaying oscillations are not observed; in the absence of enzyme input, the free enzyme concentration 
monotonically decays at an initial exponential rate. This accords with the linear analysis performed above. Though 
the form of the QSSA-based approximation suggests that the open system might oscillate in certain parameter 
regimes, multiple complementary arguments demonstrate that this is not in fact the case. 



6 Discussion 

Biochemical reactions in vivo take place under conditions in which the concentrations of substrates and active 
enzymes vary. Yet the dominant framework for modeling biochemical reaction dynamics, the Michaelis-Menten 
formalism, presumes that such reactions occur in an essentially static environment, isolated from biochemical 
inputs. This assumption is made explicit, and is indeed crucial, in classical derivations of Michaelis-Menten dy- 
namics using the QSSA [5|. Relaxing this assumption in the QSSA would therefore seem to be a natural, if 
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not essential, step in evaluating theoretically the suitability of the Michaelis-Menten formalism for modeling bio- 
chemical activity in living organisms. Prior studies have investigated the effects of both constant and time-periodic 
substrate inputs on the derivation and validity conditions for standard and reverse QSSA [ 8 , 10 1. They speak to the 
fact that in living systems, enzymatic reactions generally occur in networks where the product(s) of one reaction 
serve as substrate(s) for others, and therefore the amount of substrate available in downstream reactions is hardly 
constant as assumed in the classical QSSA. 

In the present study, we addressed the fact that in living systems, levels of enzymes in their unbound, ac- 
tivated forms fluctuate for a variety of reasons, among them changing levels of gene activation, allosteric and 
non-allosteric competition, and enzymatic degradation and denaturation. The effects of such variations are as- 
sumed to be negligible in the classical QSSA, and have remained unexamined up to now. In this work, we have 
performed the first examination, to our knowledge, of the effects of enzyme input and removal on the derivation 
and validity conditions for the QSSA. Our results further delineate the nature and suitability of the Michaelis- 
Menten formalism for in vivo biochemical reaction networks. 

We studied a particular modification of the canonical "closed" enzymatic reaction which incorporates generic 
forms of enzyme input and removal and is thus applicable for a variety of biological scenarios. Our approach 
followed the steps of the classical standard QSSA analysis, with the critical difference that we could not assume the 
total amount of enzyme (free and bound) remained constant during the reaction. The absence of this conservation 
relation distinguishes the mathematical problem considered here from all prior work on QSSA, in which the 
conservation of total enzyme is exploited in the crucial first step of reducing the dimension of the system. We used 
a different, integral form enzyme conservation relation as the starting point for our analysis of the open system, 
from which followed a recasting of the system of first order ODEs as second order differential equations. The 
reformulated system admits straightforward closed form inner solutions for enzyme and complex concentrations 
(Section [2}. These inner solutions reduce to the standard results from Q as enzyme input and removal rates go to 
zero. 

Continuing the usual steps of standard QSSA analysis, we also computed expressions for post-transient quasi- 
equilibria of enzyme and complex in terms of substrate level (Section [T2| . The expressions we obtained appear as 
intuitive corrections to the quasi-equilibrium formulas from standard Michaelis-Menten dynamics. They predict, 
in particular, the elevation of long term enzyme levels by the ratio of enzyme input to the complex formation 
rate, i.e. the equilibrium free enzyme accumulation, and the depression of long term enzyme levels by the ratio of 
enzyme removal to the complex formation rate, i.e. the equilibrium free enzyme subtraction. We note, however, 
that the corrected expressions presented here have only limited applicability due to their implicit assumption of 
static total enzyme levels, an assumption from which the open system deviates increasingly over the course of the 
reaction. Additional work, likely including novel analytical insight, will be necessary to overcome this limitation 
in order to obtain accurate outer solutions valid for the entire post- transient period of the open reaction. One 
potential avenue of investigation would be to modify the form of enzyme input in the open system ( [T0| ) such that 
it saturates, a scenario perhaps more realistically modeling the influence of varying gene expression on enzyme 
levels. 

Our examination of the pre- steady- state period was more successful: we calculated the time scale for the 
initial transient period and derived validity conditions for the inner solutions (Sections |3 . 1 1 and [33] ) . Both sets of 
expressions closely match standard QSSA results for closed reactions. The time scale for the transient period is 
shortened slightly from that of the closed system, while the primary inequality for QSSA validity is relaxed. The 
amount of both corrections depends on the rate of enzyme removal (not enzyme input), which is consonant with 
the guiding intuition behind the standard QSSA: when substrate significantly exceeds enzyme concentration, free 
enzyme rapidly binds with substrate and achieves quasi-equilibrium. Enzyme removal increases the amount by 
which substrate exceeds enzyme, and, by reducing the amount of free enzyme available, decreases the length of 
the pre- steady- state transient. 

With the derivation of inner solutions, time scale estimates, and validity conditions for QSSA for the open 
system, we thus achieved in large measure our first aim set forth in the Introduction, the extension of standard 
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QSSA to an enzymatically open reaction. We largely relied upon numerical investigations to address our second 
and third aims. Investigating the suitability of the closed system, and hence classical Michaelis-Menten formulae, 
in approximating the behavior of open biochemical reactions, we found it to be reasonably accurate so long as free 
enzyme input and removal occur several orders of magnitude more slowly than other reaction rates. While this may 
frequently be true in vivo, it is far from the rule, as described in the Introduction. Furthermore, in approximating 
the open reaction, the closed system underperforms the inner solutions we derived in (Section [2} at all rates of 
enzyme input and removal. Once the enzyme input and/or removal are at one- thousandth to one-hundredth the 
rate of other reaction steps, the behavior of the open system departs substantially from the Michaelis-Menten 
formalism. Our inner solutions accurately approximate this behavior for some time beyond the duration of the 
initial transient period, up to the point that enzyme input and/or removal occur as rapidly as the other components 
of the reaction. 

Other studies have demonstrated that the reduced systems obtained via QSSA approaches may not faithfully 
reproduce the full range of behaviors available to the original, unreduced systems (TT][T8). m particular, it has 
been shown that the canonical enzymatically-catalyzed biochemical reaction {l} admits oscillatory solutions with 
substrate input, though the QSSA-derived reduced system does not |[TT1 . In the course of investigating the range of 
dynamical behavior in the open system, we noted the apparent possibility of the opposite occurrence (Section [5]): 
though analysis of the the full open system shows that damped oscillatory solutions are impossible, the reduced 
equation |T4| ), takes the form of a forced damped harmonic oscillator and would at first glance appear to admit 
decaying oscillatory solutions. The impossibility of finding physically admissible parameter values to realize such 
solutions in either the full or reduced versions of the approximate open system, however, may be explained by 
simple arguments. 

The theoretical approach taken in this paper follows the methods and hypotheses of the standard QSSA, with 
the primary result that a more suitable inner approximation to the dynamics of enzymatically open biochemical 
reaction was derived. The reverse QSSA scenario was examined numerically; anticipated future work includes 
closer theoretical examination of reverse QSSA in open systems. An alternative mathematical approach to QSSA 
involving a different change of variables, total QSSA, has been successfully used to enlarge the range of validity 
for QSSA analysis [6, 9|. We have applied total QSSA techniques to the open system considered here, and we 
shall report our findings in another paper. 
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